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Motivation 

Advancements in long term, 
in-space, cryogenic 
propellant storage and 
transfer science and 
technologies are key to 
increasing safety, 
decreasing cost, and 
increasing payload mass of 
NASA’s space missions. 


Overall Goal 

Perform slosh experiments 
with water and LN2 to 
generate data relevant to 
benchmarking and 
expanding CFD simulation 
tools to characterize slosh 
dynamics of cryogenic 
propellants in 1g and 
microgravity storage, 
management, and transfer 
applications. 
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Progress Towards 
Goal/Agenda 


Analytical models implemented 
Ground-based test platforms 
Damping tests 

Forced sinusoidal excitation tests 

CFD simulations 

Next year plans and future work 
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Analytical Models 


• Modes 

• Wall height and forces 


d 2 V n , 2 2 D n d 2 x 

dt 2 + "" Vn ~ An “ E n dt 2 




npilaR) 2 



d 2 t] n 
dt 2 


• Assumptions: rigid sphere, inviscid, first 
three asymmetric (m=1 azimuthal wave 
number) modes 

• Implemented in MATLAB 
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Analytical Modes 
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Force Parameter vs. Excitation 
Frequency Parameter 
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TANK AND 
INSTRUMENTATION 
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Tank 


• ~30cm inner diameter 

• Aluminum 

• O-ring seals 

- Rubber for water 

- Lead wire for LN2 

• Aluminum ruler screwed to 
inside 

• Brass thermoprobe pass-thrus 

• Hole for camera 



Outside of upper half 



Inside of lower half 
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Tank 

Low density polyurethane 
foam insulation 

Aluminized nylon radiation 
barrier 

Polycarbonate mounting 
brackets 

Aluminum extrusion frame 
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Instrumentation 


PCB piezoelectric force sensor in forcing axis 

Accelerometer in forcing axis 

LVDT for position measurement 

7 Thermoprobes and 4 patch-type 
thermocouples 

National Instruments data acquisition system 

- Synchronized 

- Capable of 6kHz sampling, 1-2kHz used 
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Tank cross-section and 
thermocouple locations 



Thermoprobe 


Carbon fiber sheath 


Insulation 


Patch-type 

thermocouple 


Mounting bracket 


Patch-type 

thermocouple 


Tank flange 


Brass thermoprobe 
pass-through 
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• IDS Machine vision 
camera 

- 1 MP at 34fps 

- Fisheye lens 

- C frame grabber 

• Lighting challengi 

• Anti-distortion 
software did not 
help much 
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CFD 
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CFD Approach 


• STAR-CCM+ 

• Tank modeled as perfect sphere with a ring- 
shaped pressure outlet 

• Hexahedral dominant mesh, wall prism layer 

• Implicit unsteady, 2 nd order time and space 

• VOF formulation 

• Laminar 

• Incompressible, isothermal 

• Position-commanded motion 
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CFD Approach 


3-axis forces, 3-axis moments, and wall 
height recorded vs. time 

Wall height recorded via a field function 
that emulates the ruler inside the tank 

Tabularized position input from filtered 
LVDT data attempted 

- Ultimately used a pure sinusoidal excitation 

Computational resources primary limiter 
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Mesh size and Time step 


Dependence Study 


Meshes 


- Hexahedral: 115k, 340k, 580k, 1.3M cells 

- Polyhedral dominant: 1.2M cells 

• Time step 

- 0.001s and 0.0005s 

- 0.0002s for largest meshes 

• Test case: 50% fill, 1 .5Hz, 3mm amplitude excitation, 
10s runtime. 

• 580k mesh and 0.001 s dt selected despite mesh and 
time dependence due to computational resource 
limitations 
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RESULTS 
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Static Boil-Off Tests 


• Insulated tank, no motion 

• Filled to 90% following chill-down process 

• Fluid level measured by eye from internal ruler for 4, 
approximately 20 min periods 

• For validating a GFSSP model 

• Table starts around 86% and ends around 7% 


Arc height start [m] 

Arc height end [m] 

A Volume [m 3 ] 

Boil-off rate [kg/h] 

0.315 

0.303 

0.00054 

1.307 

0.23 

0.224 

0.000412 

0.997 

0.165 

0.157 

0.000382 

0.928 

0.132 

0.12 

0.000352 

0.851 
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Damping Test Setup 
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Damping Tests 


10 volume fractions 

2 cycle sinusoidal excitation at 1 Hz and 
2Hz, various amplitudes 

Data collected for 30-1 20s 

Data postprocessed in MATLAB 


Damping Calculations 


Force decay: 

Wave amplitude decay: 
Damping factor: 


8 , = [nil- 


n 


8wa ~ Imiz 


'Fyi + 1 
W r 


n 


8 


M 


n+ !■ 


Y = 


2n 


Peak-to-peak amplitude used 

Smoothed and averaged values 
presented 

Various correlations based on fill level 
tank radius, viscosity, and gravity 


VWg 
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First Mode Frequencies 


Water 


LN2 




h/R 
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vs. Sumner 
Experimental Data 


• 



— Eqs. (7) and (8) 


Eq. (6) 


* Eqs. (9)-(ll) 
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h/R 

• 32in diameter tank with water 

• Sense of the magnitude and variance expected 

• 400%+ difference from correlations for some points 

• XO/D for this project range from 0.0169 to 0.1014 

• Dependence on excitation amplitude apparent 
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Logarithmic Decrement 

- Water 
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Logarithmic Decrement 

- LN2 
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Eqs.(7) and (8) 


Eq.(6) 
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Eqs.(9)-(1 1) 
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Fluid surface comparisons 

-Water 
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Fluid surface comparisons 

-LN2 
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Videos 


(place holder. Too big to embed) 
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Damping Conclusions 


• Force decay method likely more accurate 

• Splashing, high nonlinearity, many other modes 
excited 

• Rotation induced in some cases 

• Clear dependence on excitation amplitude and 
frequency 

• Higher error likely a combination of: 

- Experimental error 

- Correlations not applicable 

- Correlations inadequate 

• Use better excitation for future tests 

• CFD does fair to poor job at predicting damping 
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Forced Excitation Test Setup 
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Forced Excitation Tests 

• 20% and 50% for water, 50% for LN2 

• Approximately 20 frequency/amplitude 
combinations 

- 0.5 to 4.5Hz 

• Tests repeated 3 times 

• Data collected for 30-60s depending on 
number of cycles 

• Data postprocessed in MATLAB 
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Slosh Force vs. Time 



• Water, 1.5Hz, 3.04mm, 50% 
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Slosh Force vs. Time 



• LN2, 1.5Hz, 3.04mm, 50% 
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Slosh Force vs. Time 



0 5 10 15 20 25 


Time [s] 


• LN2, 1 .59Hz, 0.93mm, 50% 
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Wave Height vs. Time 



Time [s] 


Water, 1.5Hz, 3.04mm, 50% 
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Wave Height vs. Time 



• LN2, 1.5Hz, 3.04mm, 50% 
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Wave Height vs. Time 



Time [s] 


LN2, 1.59Hz, 0.93mm, 50% 
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Force Parameter vs. 
Frequency Parameter 



• Water, 20% 
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Force Parameter vs. 
Frequency Parameter 



Water and LN2, 50% 
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Force Parameter vs. Fill Level 
and Excitation Amplitude 


0.08 r 
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Sumner65 
Water 20% 
Water 50% 
LN2 50% 

LN2 CFD 50% 
Water 50% 
LN2 50% 
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Videos 


(place holder. Too large to embed) 


Florida Institute of Technolog} 

High Tech with a Human Touch 7V 


t 


Conclusions/Results from 
forced excitation tests 


First mode resonant frequency slightly lower 
than theoretical 


Did not successfully excite the 2 nd and 3 rd 
asymmetric frequencies 

- Did seem to excite the 3 rd symmetric 

CFD does a fair job of predicting forces and 
wave height 

Rotation/swirl common at excitation 
frequencies equal to or above the first mode 
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Uncertainty/Error 


Volume measurement large source of error 

- Boiling 

- Camera resolution and angle 

- +/- 3mm to 5mm 

• At 50% fill, +/- 4mm corresponds to +/- 270mL (+- 2%) 

Tank not a perfect sphere 

Mechanical vibration noise 

No rigorous uncertainty analysis performed, 
though tests were repeated with that in mind 
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Final Conclusions and 
Future Work 


• Some confidence in the CFD models to accurately predict fluid 
slosh 

- Need to perform many more simulations 

- Larger mesh 

• Need to perform an uncertainty analysis 

• Hardware improvements necessary 

- Reading fluid heights from 100000’s of images is not feasible 

- Lighting 

- Thermocouple instrumentation 

• Damping tests need to be rerun with an emphasis on only 
strongly exciting the first mode 

• Rotational modes are high amplitude/low decay 

- Clearly important to understand, but nothing planned 


Florida Institute of Technology 

High Tech with a Human Touch 7M 


Final Conclusions and 
Future Work 


Forced excitation tests: additional volume fractions 
Fluid management devices 

- Baffles 

Free pitching axis 

“Floating tank” approach versatile 

- Free translation tests 
“In-space” part 

- Parabolic aircraft flight experiments 

- CFD modeling of past work 

• Drop tower 

• FIT-MIT SPHERES SLOSH ISS experiment 

Brainstorming possible fluid transfer experiments 
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EXTRA 
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Thermocouple Example Plot 
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Interesting Videos 


(too large to embed) 
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Mesh Cross-sections 
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Cool pictures of seal test 





